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ABSTRACT 


! 

|l 

The  essential  features  of  zonally-averaged  energy  budget  models  of  the 
earth's  climate  are  reviewed  and  abstracted  to  give  an  "effective  potential" 
for  the  appropriate  climate  variables.  The  effective  potential  permits  the 
multiple  valued  insolation  versus  ice  latitude  behavior  common  in  such 
models  to  be  considered  from  the  point  of  view  of  a smoothly  varying  energy 
function.  Next  we  view  weather  as  a short  time  scale  driving  mechanism  for 
the  long  time  scale  variations  of  climate  and  give  a formalism  for  explicitly 
realizing  a climate  effective  potential  by  integrating  out  or  averaging 
over  the  short  time  scale  weather  variations.  The  effective  potential  so 
realized  embodies  many  of  the  features  abstracted  from  the  energy  budget 
models. 
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I INTRODUCTION 


In  the  past  several  years  there  has  been  an  active  consideration  of 
models  of  the  growth  and  recession  of  the  earth's  ice  caps.  These  models 
have  their  bases  in  equilibrium  states  dictated  by  an  energy  balance 
between  insolation,  radiation  from  the  earth  after  differential  absorption 
by  ice  covered  or  by  ice  free  regions,  and  energy  transported  poleward  by 
the  atmosphere  and  oceans.  The  pioneering  work  in  these  "energy  budget" 
models  was  done  by  Budyko  (1969)  and  Sellers  (1969).  It  has  been  followed 
by  a number  of  papers  including  Schneider  and  Gal-Chen  (1973),  Gal-Chen 
and  Schneider  (1976),  North  (1975a,  1975b),  Held  and  Suarez  (1974),  Sellers 
(1973),  Su  and  Hseih  (1976),  and  Lindzen  and  Farrell  (1977). 

The  continuing  interest  in  these  simplified  models  stems  from  the  fact 
that  they  consistently  exhibit  a great  sensitivity  to  small  changes  in  the 
solar  constant.  Precisely  what  level  of  sensitivity  is  to  be  believed 
remains  a matter  of  active  discussion,  Sellers  (1977)  and  Dickinson  (1977); 
nevertheless,  since  it  is  common  for  these  models  to  permit  an  uncontrolled 
growth  of  the  earth's  ice  caps  after  a few  percent  decrease  in  the  solar 
constant,  the  significance  of  pursuing  them  is  clear.  The  physical  mechan- 
ism for  this  instability  is  not  hard  to  identify:  as  one  lowers  the  solar 
constant  the  average  temperature  at  any  latitude  decreases  and  the  ice  pack 
grows.  The  ice  pack  has  a higher  reflectivity  (albedo)  than  the  sea  or 
land  so  more  sunlight  is  reflected  which  In  turn  tends  to  cool  off  the 
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earth,  form  more  ice,  etc.  This  instability  is  countered  by  the  transport 
of  energy  poleward  by  the  atmosphere  and  oceans.  When  that  transport  is 
insufficient,  the  ice  cap  will  grow.  A review  of  this  mechanism  in  the 
general  context  of  climate  models  is  given  by  Schneider  and  Dickinson  (1974). 

In  this  paper,  after  a review  of  the  essential  points  of  energy  budget 
models,  I discuss  how  one  might  interpret  these  models  within  the  general 
framework  of  theories  such  as  those  discussed  by  Hasselman  (1976).  In  such 
theories  the  climate-weather  system  is  divided  into  short  time  scale  pheno- 
mena called  weather  and  long  time  scale  phenomena  called  climate.  The 
time  scales  involved  in  weather  may  be  as  short  as  a few  days  or  weeks; 
while  climate  time  scales  range  from  several  months  or  years  to  many  thou- 
sands of  years  in  the  case  of  ice  ages.  It  is  not  absolutely  certain  that 
one  can  cleanly  divide  the  time  scales  of  weather  and  climate  since  there 
may  be  a very  large  number  of  important  time  scales  ranging  from  days  to 
thousands  of  years;  though  it  does  seem  to  be  generally  agreed  upon, 

Lorenz  (1977),  that  such  a division  is  attractive. 

Our  procedure  here  will  take  the  energy-budget  models  and  cast  them 
into  the  language  of  an  "effective  potential"  whose  smooth  behavior  as  a 
function  of  external  parameters,  such  as  Insolation,  can  give  rise  to 
apparently  discontinuous  or  multiple  valued  equilibrium  states.  For  the 
ice  cap  problem  the  potential  will  be  a function  of  the  cosine  of  the  ice 
cap  latitude,  y,  and  the  Insolation,  q,  normalized  to  current  values.  The 
potential  V(q,y)  will  be  represented  by  a polynomial  in  y with  coefficients 
dependent  on  q (and  thus,  If  q depends  on  time,  dependent  on  time): 


2 


V(q,y) 


+ 


+ 


+ 


(1) 


« a.(q)y2 

2 


b(q)y3 

3 


c(g)y4 

4 


Three  terms  will  suffice  for  the  problem  at  hand.  The  time  dependence  of 
y will  come  from 


Iz  « -avj&ty) 

3t  3y 


(2) 


which  accounts  for  the  name  "effective  potential".  Clearly  the  inflection 
points  of  V(y)  with  positive  curvature  as  well  as  the  end  points,  y = 0 or  1, 
with  negative  slope  are  the  possible  equilibrium  states.  As  q varies,  the 
coefficients  a,b,c,  . . . will  vary  smoothly  changing  the  character  of  the 
equilibrium  states  in  a dramatic  way. 


The  nature  of  V(q,y)  will  come  from  a study  of  the  equilibrium  states 
a'la  Bud yko- Sellers.  Once  we  have  V,  however,  the  time  dependence  of  y(t), 
that  is,  the  climate,  is  rather  easily  extracted  from  the  elementary 
equation  (2).  The  potential  itself  is  an  abstract  formulation  of  the 
energy-budget  models;  its  usefulness  is  then  in  the  study  of  the  dynamic 
properties  of  climate.  In  the  ice  cap  example  the  effect  of  time  variations 
in  q on  y(t)  are  simply  followed  by  solving  (2).  So  if  one  imagines 
variations  in  insolation  to  be  at  the  source  of  ice  ages,  various  formula- 
tions of  this  hypothesis  may  be  straightforwardly  examined. 


In  addition  to  providing  an  effective  potential  I will  give  a discus- 
sion of  its  derivation  from  the  interaction  of  short  and  long  time  scale 
phenomena.  Essentially  I will  treat  a two  time  scale  system  where  weather 
is  represented  by  a gausslan  statistical  field.  Averaging  the  climate 
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(long  time  scale),  which  is  driven  by  the  weather,  over  the  weather  variables 
(that  is,  "integrating  out"  the  weather)  yields  an  effective  potential  for 
climate  with  features  similar  to  those  abstracted  from  energy  budget  models. 

• 

» 

I 

> 

> 
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II  ENERGY  BUDGET  MODELS  FOR  ICE  CAPS 


In  this  section  will  be  a quick  review  of  energy  budget  models  of 

Budyko-Sellers  type.  The  treatment  follows  North  (1975a)  in  spirit.  We 

consider  a zonally  averaged  earth  with  a surface  temperature  T depending 

on  latitude  0 and  time  t . The  earth  has  ice  caps  extending  to  6^  and 

open  sea  or  land  lying  between  0q  and  -0  . The  ice  is  assigned  an  albedo 

* 0.6  for  | © | > 6 q , and  the  ice-free  earth  is  assigned  a lower  albedo 

a * 0.4  for  | © | < © . The  net  incoming  radiation  absorbed  will  then  be 
o 


A(0 , 0 o [ t ] ) = QS(0) 


1 - a,  0 > 0 
i o 

1 - a 0 < 0 


where  we  consider  only  0 l 0 < ir/2  because  of  symmetry  about  the  equator. 
In  (3)  we  have 


Q * global  average  insolation  - 1/4  solar  constant  - 0.35  kwatt/in  , 


Qs(0)  * annual  mean  insolation/unit  area,  so  i d0  cos  0 s(0)  - 1 . (4) 


After  absorption  this  radiation  is  either  reradiated  or  transported 

poleward  by  the  ocean  and  atmosphere.  All  workers  seem  to  agree  that  the 

4 

real  black  body  radiation,  oT  , is  very  well  represented  by  the  linear  form 

4 

A + BT,  which  is  essentially  an  expansion  of  aT  about  T - 0°C  and  is 

acceptable  since  changes  in  T are  not  large.  The  coefficients  are 

2 2 
A * 200  watts/m  and  B s 1.5  watts/m  -°C. 


The  representation  of  the  heat  transport  by  ocean  currents  and  by  the 
atmosphere  is  complicated  in  detail.  It  seems  quite  reasonable,  however, 
to  adopt  the  usual  form  for  diffusion 


cos 


Je  (cos0  D H) 


(5) 


where  D is  a diffusion  coefficient  sometimes  taken  to  depend  on  6. 


All  together  these  considerations  lead  to  the  diffusion  equation 


3T(9,t) 


at 


energy  absorbed  - energy  radiated  + energy  transported 


= Qs(0)A(0,eo[t])  - (A  + BT[6  ,t])  + Jq  ( cos0  D T[0,t]) 


(6) 


where  C is  some  "heat  capacity"  and  really  serves  to  set  the  time  &._ale 
on  which  the  system  approaches  equilibrium. 

Taking  D to  be  a constant,  North  (1975a)  is  able  to  provide  an  anal- 
ytic solution  to  the  time  independent  version  of  (6).  That  is  to  say,  he 
seeks  the  possible  equilibrium  states  governed  by  the  energy  balance  indi- 
cated on  the  right  hand  side  of  (6).  The  boundary  conditions  used  are 


1,  (l-K2)* 

O X 


x = sin  0 


x=0 

and 

x=l 


that  is,  no  heat  is  transported  across  the  equator  or  poles. 


2)  lim  T(x  +e)  “ lim  T(x  -e)  = T , x ■ sin  0 


e+0 


e-*-0 


that  is,  the  temperature  at  the  ice  front  (at  9 or  x ) is  the  same 

o o 
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approached  from  below  or  above  in  latitude.  This  temperature  is  taken  to 
be  a convenient  number,  say  -10°C,  where  water  is  likely  to  be  frozen. 


» S 


x=x  + 0 
o 


3T 

3x 


x»x  - 0 
o 


that  is,  no  net  heat  is  transported  across  the  ice  front. 


The  result  of  this  kind  of  model  calculation  is  quite  clearly  presented 

in  terms  of  a graph  of  x = sin  0 versus  q , the  solar  constant  normalized 

o o 

to  the  contemporary  value  (see  Figure  1).  This  figure  is  sketched  from 
North  (1975a)  but  is  similar  to  that  found,  for  example,  by  Held  and  Suarez 
(1974).  As  demonstrated  explicitly  by  North  (1975a)  and  by  Su  and  Hsieh 
(1976)  the  regions  with  dq/dxQ  > 0 are  stable  in  the  sense  that  small 
changes  in  q lead  to  exponentially  damped  changes  in  x^(t)  computed  from  (6). 
When  dq/dxo  < 0,  (6)  has  solutions  for  xQ(t)  which  grow  until  the  ice 
front  reaches  the  equator,  xq  = 0. 


, These  models  are  quite  sensitive  to  the  connection  albedo  and  temper- 

ature or  equivalently  the  magnitudes  of  ct^  and  a.  Clearly  the  larger 
becomes,  the  more  reflection  there  will  be  from  ice  and  the  stronger  will 
I be  the  driving  mechanism  for  unstable  ice  front  growth. 

There  are  many  serious  simplifications  in  this  kind  of  model.  Perhaps 
I the  most  blatant  are  the  treatment  of  the  atmosphere  as  homogeneous  and 

essentially  cloud  free  and  the  use  of  quantities  averaged  over  longitude  as 

» 

well  as  over  seasons.  Corrections  such  as  these  seriously  complicate  the 
I calculations  in  the  model  and  yet  result  in  rather  similar  patterns  for  the 

i; 


» 
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Figure  1 Coetne  or  tine  of  the  ioe  from  latitude  m a function  of  itoiatton.  Thh  curva  la  mora  or  taa  takan 
from  North  (1975  a),  but  it  is  rapraaantathra  of  energy  budget  modai  caiculations. 


I 


I 

t 

I 

’ 

! 


» 
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q versus  xq  behavior,  Sellers  (1977).  For  the  present  we  will  accept  the 
results  of  Figure  1 as  representative,  if  perhaps  absent  of  fine  structure, 
and  attempt  to  give  a general  framework  for  these  models.  This  occurs  next. 


» 


I 


I 
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III  GENERALITIES  ABOUT  THE  ENERGY  BUDGET  MODELS  : EFFECTIVE  POTENTIALS 

We  begin  by  examining  Figure  1.  When  q is  very  large,  q > q^  * 1.3, 

there  is  a single  allowed  equilibrium  state  which  has  y = cos  0 =0,  that 

o o 

is,  no  ice  cover  at  all.  When  q is  lowered  to  below  q^  but  above 

q2  * 1.10,  we  have  three  possible  equilibrium  states:  y^  = 1 (ice  covered 

earth),  * 0,  and  a set  of  intermediate  values  depending  on  q . The 

states  yQ  = 0 or  1 are  stable;  the  other  is  unstable.  When  q goes  below 

q0  * 1.10  but  is  above  q_  * 0.94,  there  are  still  two  stable  states  y = 1 

^ j o 

and  yQ  in  the  region  0 to  0.8.  The  present  ice  cap  has  yQ  * 0.32  and 

clearly  lies  in  this  region.  There  is  also  a possible  unstable  state  in 

this  regime.  Finally,  when  q < q^  » 0.94,  there  is  only  one  possible  equil- 
ibrium state  y ■ 1,  and  the  ice  cap  grows  in  an  uncontrolled  way  to  cover 
o 

the  earth. 


I 

l 

i 


1 

< 


» 


» 


* 


Such  a situation  is  reminiscent  of  a system  undergoing  various  phase 
transitions.  The  external  parameter  in  most  phase  tranp  '.tions  is  the  temp- 
erature, or  the  concentration  of  some  impurity  or  perhaps  the  magnitude  of 
an  external  magnetic  field.  Here  we  have  the  possible  phases  of  the  system 
governed  by  the  insolation  q . In  the  theory  of  phase  transitions  it  has 
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proved  convenient  to  characterize  the  system  by  an  effective  free  energy 
taken  to  be  a polynomial  in  a mean  field  and  whose  coefficients  depend  on 
the  external  parameter.  Such  a phenomenological  approach  was  initiated  by 
Landau  and  co-workers  (see  Landau  and  Lifshitz,  1968)  and  has  proven 
especially  useful  as  a compact  description  of  an  intrinsically  complicated 
system,  for  example,  a ferromagnet. 


In  the  present  case  we  want  to  think  of  yQ(t)  as  a "mean  field".  In  a 
more  realistic  model  the  ice  front  would  depend  on  the  longitude  X and 


the  mean  field  would  be  yQ(X,t),  as  computed  from  T(0,X,t)  = Tq.  We  now 
assign  an  effective  potential  V(q,yQ)  to  the  system  and  determine  the  time 
dependence  of  yQ(t)  from 


dyQ(t)  3V(q,yQ) 


In  the  present  case  it  will  be  sufficient  to  choose 


v(q.yQ)  = \ a<9> y20  + J b(q)y^  + \ c(q)y£ 


"dt^  = _y0(t:)[a(q)  + b(q)yQ  + cCq)y^| 


To  describe  the  situation  in  Figure  1 we  need  the  following 


1)  When  q > q^ 


a(q),  b(q) , c(q)  > 0 . 


This  means  V(q,yQ)  appears  as  in  Figure  2.  Whatever  the  initial  value  of 

y it  will  "roll  down  the  hill"  toward  y ■ 0 (no  ice)  at  a rate 
o o 


• —v 


V 

exp  -/  a(q[t'])dt'  for  large  t . 

2)  When  q reaches  some  value  qQ  > q^,  the  coefficient  c(q)  becomes 
negative  and  at  q « q1  the  slope  of  V(q,yQ)  at  yQ  - 1 is  zero.  For  q s.  q.^, 
but  q > q2  we  have 


3V 

a(q),  b (q)  > 0 , c(q)  < 0 , and  -r— - 

3yo 


< 0 , (11) 


yo  * 1 


as  shown  in  Figure  3.  If  y > y for  q_  < q < q, , the  ice  cap  will  grow  in 

o i.  X 

time  to  y =1.  If  y < y,  the  ice  will  vanish  in  time.  The  point  y * y 
o o o 

is  a possible  equilibrium  state,  though  it's  clearly  unstable  since 


32V 


3y 


< 0 


o o 


3)  When  q < there  develops  a stable  equilibrium  yQ  away  from 
yQ  *=  0.  In  the  domain  q^  < q < q2  we  have 


a (q)  , c (q)  < 0 , and  b(q)  > 4a(q)c(q) 


(12) 


so  the  roots  of  3V/3yQ  , 


-2c(q)y+ * b(q)  +.  /b(q)  - 4a(q)c(q) 


are  real  and  distinct.  Now  the  potential  appears  as  in  Figure  4. 


If  the  starting  value  of  yQ  is  less  than  y+  , then  yQ(t)  will  move 
toward  yQ  m y_  as  time  goes  on.  The  rate  of  approach  will  be  governed  for 
long  times  by 


32V 


3y 


o 'o 
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The  value  y * y+  Is  a possible  equilibrium  value,  but  again  from  Figure  4 
is  clearly  unstable. 

4)  Finally  when  q < , the  transport  of  heat  poleward  is  not 

sufficient  to  overcome  the  ice  albedo  instability  and  the  ice  grows  to 

y =1.  In  this  case 
o 

a(q) , c (q)  < 0 and  b(q)2  < 4a(q)c(q)  . (15) 

Figure  5 shows  this  growth  of  y^  , whatever  its  initial  value,  in  graphical 
form. 


A possible  simple  parametrization  of  a,  b,  and  c is 

a(q)  = a(q-q2) 

b (q)  «=ye(q-q3)  + 4ay  (q-q2)  (q-qQ)1  (16) 

c(q)  » y(q-q0) 

with  a,  8,  y > 0 , qQ  > q^  and  q^  defined  by 

- 0 - a(q3)  + b(q3)  + cfq^  (17) 

V1 

q-qx 

The  introduction  of  the  effective  potential  V(q,y  ) has  thus  far  only 

o 

been  a transcription  of  Figure  1 into  other  terms.  It  proves  useful  when 
we  wish  to  follow  the  time  dependence  of  yQ  by  solving  (9) . Taking  q to 
be  time  dependent  we  need  to  solve 
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(18) 


dy  (t)  , , 

dt “ -yo(t)  + Mt)yQ(t)  + c(t)yQ(t)  J 

for  the  parametrization  of  a,  b,  and  c in  (16)  and  for  the  chosen  time 
dependence  of  q . It  is  clear  that  as  q changes  we  move  from  regime  to 
regime  in  V(q,yQ)  and  even  small  changes  in  q can  result  in  significantly 
different  equilibrium  states.  We  appear  today  to  be  in  the  regime  (see  Fig.  4) 
q > q^,  y < y+  ; if  q were  to  dip  below  q^,  the  ice  caps  would  grow  from 
their  present  value.  If,  during  this  growth,  q moved  back  to  be  > q^, 
but  at  that  time  y^  were  greater  than  y+,  the  ice  caps  would  keep  on 
growing.  If  y were  less  than  y+  when  the  switch  occurred,  then  y would 
go  back  to  y . Following  yQ(t)  for  any  choice  of  q(t)  involves  the  rather 
straightforward  solution  of  (18) . 

This  equation  is  known  as  the  Abel  equation  of  the  first  kind  (Murphy, 
1960,  p.  23-25)  and  is,  in  general,  not  soluble  in  closed  form.  However, 
near  b 1 0 or  c * 0 it  may  be  easily  solved.  As  an  example  of  an  oversim- 
plified ice  cap  model  consider  a system  where  a(q)  = a^q(t)  - q^j,  b(q)  = 0, 
and  c(q)  * y ^q(t)  - q£j  with  9a  > 9C  • When  q > the  system  has  a single 
stable  state  at  yQ  = 0 ; when  q < qc,  y^  = 1 is  the  only  stable  state.  For 
qc  < q < q^,  various  Intermediate  states  exist.  In  this  situation  yQ(t)  is 
easily  found  from  (18)  to  be 


yQ(0)  exp  - J a(x)dx 


y0(t) 


1 + 2y‘(0) 


t 

I F(x)dx 
o 


(19) 
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F(x)  = ^exp  - 2 J du  a(u^  c(x) 


When  c > 0,  a > 0,  yQ(t)  tends  to  yQ  = 0 for  t large.  When 


c > 0,  a < 0,  y tends  to  s 
o 


. When  a and  c are  both  negative. 


yQ  grows  until  yQ  = 1.  This  can  all  be  qualitatively  traced  by  drawing  the 
effective  potential  for  these  cases.  The  implication  of  any  hypothesis  for 
the  time  dependence  of  the  insolation,  q(t),  on  the  growth  and  recession  of 
the  ice  caps  can  now  be  directly  drawn  from  (19)  and  (20) . 


» 


I 


> 


The  use  of  the  effective  potential  is  clearly  not  limited  to  the  over- 
simplified zonally  averaged  models  considered  to  this  point.  Whenever  we 
have  a climate  variable  (J>(x,t),  we  may  find  a potential  for  it  from  the 
underlying  dynamics  or  from  a study  of  the  equilibrium  states  where  = 0 
as  in  the  Budyko- Sellers  approach.  The  time  dependence  will  then  be  given 
by 

3d)  “3V 

— (x,t)  * — (4>,  external  parameters)  , (21) 

and  the  various  equilibrium  states  will  be  classified  by  the  solutions  of 

2 2 

dV/d<t>  m 0 and  their  stability  determined  by  the  sign  of  3 V/3<t>  at  these 
inflection  points. 


IV  WEATHER  AS  A FEEDBACK  MECHANISM  FOR  CLIMATE 


In  this  section  we  give  a formalism  which  allows  one  to  determine  the 
behavior  of  long  time  scale  variables  (climate)  in  interaction  with  a set 
of  variables  whose  time  scales  are  much  shorter  (weather).  This  idea  has 
been  stressed  by  Hasselmann  (1976)  who  has  discussed  a "slow"  system  in  the 
presence  of  a "fast"  statistical  background  with  no  feedback.  Feedback,  or 
non-linear  effects,  are  essential  to  the  Issues  considered  above,  so  we  will 
have  to  treat  them  here. 


The  idea  is  that  there  is  a field,  W(x,t),  of  weather  variables  whose 
time  scale  and  perhaps  length  scale  are  small  compared  to  those  of  another 
field  <(>(x,t)  of  climate  variables.  Many  variations  of  W(x,t)  will  occur 
during  any  significant  motion  of  <(>(x,t),  and  we  will  assume  that  W(x,t) 
acts  as  a statistical  background  for  the  variation  of  4>(x,t)  and  further 
that  W(x,t)  is  a random  variable  at  each  point  of  space,  x,  and  time,  t, 
with  a gaussian  distribution.  The  field  <t>(x,t)  will  be  driven  by  its  inter- 
action with  W(x,t),  and  our  goal  is  to  find  the  behavior  of  4>(x,t)  once  the 
variations  in  W(x,t)  have  been  averaged  out. 

To  begin  we  require  a description  of  the  weather;  our  rapidly  varying 
phenomenon.  We  give  this  in  terms  of  a Lagranglan  density  for  the  field 
W(x,t) 
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L 

w 


2 3W(x, t)  3W(x,t)  _ t_2  ~ + 3W(x,t)  9W(x,t)  + 

1 at  at  2 [ u,c;  at  at  wu,t; 


VW (x , t ) 


VW(x,t) 


a W(x,t)  W(x,t) 

W 


(22) 


In  this  Lagrangian  the  coefficients  x^  and  are  time  scales  for  the  oscil- 
lation and  diffusion  of  weather,  and  1 is  a length  scale  governing  weather; 

w 

a^  is  a dimensionless  number.  All  of  these  coefficients  are  in  general 
functions  of  external  parameters  such  as  the  insolation.  In  the  Lagrangian 
appears  not  only  the  weather  field  W(x,t)  but  also  a mirror  image  weather 
W(x,t)  which  is  necessary  because  of  the  dissipative  term  proportional  to 
x 2 (Morse  and  Feshbach,  1953,  Section  3.3).  The  dynamics  of  weather  come 
from  the  extremum  of  the  action 


S (W,W)  = (dD  x dt  L 

w J w 


(23) 


where  D is  the  number  of  space  dimensions.  Variations  of  yield 


t2  + x , j-  - £2  V2  + a ) W ( x , t ) = 0 
V 1 3t2  2 8t  w w/ 


(24) 


and 


-^2  ~ T " *2  V2  + a ) W ( x , t ) = 0 
\ 1 3t2  2 8t  w w/ 


(25) 


which  exhibits  the  different  roles  played  by  W and  W. 


Only  quadratic  terms  appear  in  Lw  since  we  are  assuming  that  weather 
acts  as  a gaussian  random  field  (see  below) . In  some  sense  the  non- 
linearities  of  weather  have  been  hidden  in  the  term  proportional  to  x^  in 
since  this  allows  oscillations  of  the  system.  If  one  is  truly  interested 
In  the  weather,  then  additional  non-linearities  may  be  represented  by 
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C 


quartic  and  higher  order  terms  in  L 


We  want  to  treat  W(x,t)  as  a statistical  variable  at  each  point  of 
space  and  time.  To  this  end  we  must  assign  a probability  density  to  this 
distribution.  It  is  natural  to  choose  this  to  be 

P (w(x,t),  W ( x , t ) ) = exp  - S (W,W)  (26) 

so  that  the  most  probable  situation  is  that  which  has  minimum  action;  that 
is,  the  state  of  the  system  which  follows  Newton's  laws  is  most  likely. 
Furthermore,  with  only  quadratic  terms  in  L^,  this  probability  density  now 
describes  a set  of  gaussian  random  variables.  Averages  of  any  function 
F(W,W)  are  given  as  usual  as 


<F(W,W)> 


n 


JdW(x,t)  dW(x,t)  F (W,W) 
x!t  /dW(*»t)  dW(x,t)  Pw 


P 

w 


For  example,  the  "average  weather"  has  been  chosen  so 


(27) 


<W(x,t)>  = <W(x,t)>  = 0 

So  if  the  weather  variable  under  consideration  were  the  temperature  T(x,t), 
we  would  be  setting  W(x,t)  ■ T(x,t)  - T , with  T the  actual  mean  temperature. 


The  correlation  function  for  weather  variables  is,  however, 

<W(x1,t1)  W(x2,t2)>  - <£(5^,^)  W(x2,t2)>  - 0 , (28) 

and 

<W(x1,t1)  W(x2,t2)>  - | Cw  (Xj-Xj,  ^-tj)  (29) 
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with 


(T1  h + t2  FT  - K + - - S°W6M 

0 t 


(30) 


so  C (x,t)  is  the  Green  function  for  the  weather.  The  derivation  of  these 
w 

results  follows  from  the  multidimensional  gaussian  integral 


/a  dvj  -[1].1viA«vJ  + AVj] 

■ irn/2  (det  A)"1/2  exp  j l B (A-1)  B 

4 i,j=l  1J  J 


(31) 


with  the  labels  i and  j being  continuous  as  is  appropriate  for  fields 
v = v(x,t).  Averages  are  conveniently  evaluated  by  taking  derivatives  with 
respect  to  B . 


This  weather  will  drive  our  climate  variables,  so  we  want  to  describe 
the  climate  first  in  the  absence  of  short  term  variations  and  then  in  inter- 
action with  them.  Without  weather  we  take  for  a climate  Lagrangian  density 

L,  = Ti  $(*»t)  <Kx,t)  - -y  ^(x,t)  (x,t)  - ||-  (x,t)<t>(x,t)J 

- A2  7|(x,t)  • V<J>(x,t)  - Ac  i(x,t)if(x,t)  - y $ (x,t)2<Kx,t)2  , (32) 

and  for  the  interaction  between  climate  and  weather  we'll  take 
L+  “ - [W(x,t)  + W(x,t)]  [4>(x,t)  + $(x,t)] 

- *2  [W(x,t)  + W(x,t)  ] $ (x,t)4>  (x,t) 

- X3  W(x,t)  W(x,t)  [$(x,t)  + $(x,t)]  . (33) 
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The  assumption  that  weather  is  gaussian  restricts  us  to  quadratic  terms  in 
W and  W ; higher  order  polynomials  in  $ and  | may  appear.  Symmetry 
between  W and  W and  between  and  $ is  maintained  consistent  with 

their  interpretation  as  diffusing  and  anti-diffusing  fields;  see  Equations 
(24)  and  (25). 


At  this  stage  we  have  a set  of  coupled  fields  where  clearly  the  dyna- 
mics of  the  rapidly  varying  W drives  <J>  as  we  see  in  the  Euler-Lagrange 
equations  for  <f> 


(Tl77+T2H-'c,2  + Ac)»S't)  ’ 

-[X^  + X2  $(x,t)]  [W(x,  t)  + W ( x , t ) ] - X3  W(x,t)W(x,t)  - Ai(x,t)2<Kx,t)  . (34) 


The  usual  Brownian  motion  treatment  of  climate  as  in  Hasselmann  (1976)  (see 
Wang  and  Uhlenbeck,  1945),  sets  X^  a X^  = 0 and  using  (34)  derives  dynamic 
equations  for  <$(x,t£>  or  ^(x.t)  )>  by  utilizing  the  statistical  properties 
of  W . The  interaction  or  feedback  represented  by  the  X2  and  X^  terms  are 
important.  In  the  ice  front-insolation  questions  of  before  they  are  the 
effective  transport  or  stabilization  mechanisms  for  the  climate  <J>  . The 
term  in  X^,  for  example,  provides  a space-time  dependent  contribution  to 
the  balance  between  the  oscillation  and  diffusion  terms  and  directly  affects 
the  stability  of  climate  motion. 


Through  its  Interaction  with  statistical  weather  the  climate  becomes 
a statistical  quantity  Itself.  The  statistical  weight  for  the  interacting 
climate-weather  system  will  be 
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p 


P(W,W,$,$)  - exp  - S(W,W, <{>,$) 


(35) 


where  S is  the  full  action 


S(W,W,*,$)  = fd°x  dt[L  + L + L ] 

J w c lJ 


(36) 


The  statistical  properties  of  the  climate  alone  are  given  by  the  effective 
action 


exp  - Sgff  (<M) 


J x,t 


dW(x, t)dW(x, t)  exp  - S(W,W,<t>  ,$)  (37) 


which  amounts  to  integrating  out  the  weather  variables.  The  statistics  of  $ 
are  now  in 


with 


where 


Seff  = /d°x  dt  Leff 

Leff  = Lc  (4>  »$ ) - y fdPydr  U^(x,t),  <j>(x,t)) 

* Gw(*-y.  t-T)  u(<Ky,x),  $(y,T)) 

UW,^)  = X^iKx.t)  + $ (x,t)]  + X 2^  (x, t)$  (x,  t)  , 


(38) 


(39) 


(40) 


and  GW  is  the  interacting  weather  Green  function  satisfying 
2 

\ f~2  + t2  h ‘ ll  v2  + aw  + X3  (^(x*t)  + *<*»*>) 

9 1 


V*’t) 


- «u(x)«(t) 


(41) 


Looking  at  (39)  we  see  that  the  weather  has  provided  non-linear  feed- 
back terms  for  the  climate  variation  which  are  apparent  in  the  terms  of 
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L ,,  with  no  space  or  time  derivations,  which  we  call  the  effective  potential 
et  t 


V2 


V ,,($,  <t>)  = A H + (4>  + *)G  (<t>  + *) 

err  c 2 w 

[(4.  + $)Gw$$  + MGw(4  + $)]  + |$4>(A  + , (42) 


using  a more  or  less  obvious  shorthand  notation. 

The  effective  potential  is  the  key  result  of  this  section.  It  shows 

how  the  weather,  when  "removed"  from  the  climate  problem  by  integration 

over  dWdW,  has  provided  non-linear  feedback  terms  in  the  potential  governing 

the  motion  of  <p  > the  climate.  To  exhibit  the  connection  with  the  previous 

section  a bit  more  explicitly,  let's  look  more  closely  at  G^.  The  space 

and  time  scales  in  the  integration  required  for  (39)  or  (42)  are  climate 

space  and  time  scales  and  are  larger  than  or  or  If  we  simplify 

a bit  more  by  neglecting  X ^ (<p  + 4)  in  (41)  as  well  as  the  space  and  time 

derivatives  because  of  the  previous  argument,  we  may  approximate 

G (x  - y,  t - t)  by 
w 


G (x  - y,  t - t)  = - ~~  6°(x  - y)6 (t  - t) 

w a 

w 


(43) 


This  formula  is  accurate  as  long  as  a i4  0.  Next  imagine  the  climate  dissi- 


w 


pation  term  in  is  small  so  4 ■ ij>.  The  effective  potential  becomes 

v(<i,*>  - ^ *2  + ^ «3  + -c-^-  ** 


(44) 


with 


iM 

2 


2A.(qr 

V>  - Tw 


(45) 
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(46) 
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c(q) 

4 


x1(q)x2(q) 

a (q) 
w 

1 / X2(q)i 

5 AW  -77B 


(47) 


where  we  have  restored  the  dependence  of  the  coefficients  on  the  external 
parameter  q . 


The  discussion  of  the  previous  section  can  now  be  carried  out  as 

before.  The  climate  variable  will  be  identified  with  y^,  and  the  dynamics 

of  yQ  are  as  described  above.  The  value  of  the  considerations  here  are 

primarily  to  give  a framework  for  the  effective  potential  formulation.  One 

is  able  to  identify  directly  the  parameters  a(q),  etc.  of  the  effective 

potential  in  terms  of  the  parameters  governing  the  underlying  dynamics  of 

short  term  phenomena  which  drive  the  climate.  Additionally  a certain  insight 

is  provided  for  the  meaning  of  the  effective  potential.  The  dynamics  of  the 

climate  are  to  be  viewed  as  a most  probable  path  in  <t>  space;  namely  that 

path  which  minimizes  the  effective  action  S The  full  motion  is  much 

ef  f 

more  complicated  than  given  by  (44)  since  it  involves  the  high  order  non- 
linearities  contained  in  G^.  One  may  proceed,  however,  with  some  hope 
that  the  approximations  to  will  suffice. 


: 
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V DISCUSSION 
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In  this  note  we  have  fitst  of  all  tried  to  isolate  some  general 
features  of  energy  budget  models  of  the  Budyko-Sellers  variety  by  casting 
them  into  an  effective  potential  or  effective  free  energy.  Extrema  of  this 
free  energy  with  respect  to  the  climate  variables,  <f>  , will  yield  the 
equilibrium  states  of  the  energy-budget  models,  and,  as  explained,  above 
the  time  dependence  of  the  models  may  be  directly  read  off  from  the  equation 

~ = - (Effective  Potential) 

0 1 0$ 

In  the  ice  front  problem  tf>  is  taken  to  be  the  cosine  of  the  ice  front 
latitude,  yQ(t),  and  the  dynamics  then  lie  in  a first  order,  non-linear 
equation  yQ(t).  Such  equations  may  sometimes  admit  of  analytic  solutions 
and  may  be  easily  treated  numerically.  Indeed  the  major  use  of  the  effec- 
tive potential  viewpoint  is  precisely  the  ease  with  which  the  time  dependent 
problem  may  be  studied.  In  addition  it  allows  one  to  see  developing  in  a 
smoothly  varying  manner  the  apparent  discontinuous  behavior  in  the  y^  versus 
q plots  typified  by  Figure  1.  As  noted  in  the  text  this  is  just  the 
situation  in  the  theory  of  phase  transitions  where  a smoothly  varying  free 
energy  may  give  rise  to  an  apparently  discontinuous  phenomenon  such  as 
freezing  of  a liquid  or  spontaneous  magnetization  of  a ferromagnet. 

The  second  part  of  our  discussion  has  concerned  the  framework  In  which 
to  place  the  effective  potential.  Treating  weather  as  a short  time  scale 
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gausslan  random  field,  we  were  able  to  construct  an  effective  action  for 
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> 
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climate  (long  time  scale  phenomena)  in  interaction  with  weather;  i.e., 
driven  by  it.  The  key  point  was  to  assign  a probability  density  to  the 
weather  field 


P = exp  - (classical  action)  , 

where 

classical  action  = yd(space)dt  [Langrangian  density  for  the  field]  . 

In  analogy  with  usual  statistical  physics  we  would  then  interpret  the 
action  as  an  "entropy"  equal  to  - log  (Probability)  . The  most  probable 
state  extremizes  the  action  and  is  governed  by  the  Euler-Lagrangian  equations 
(Newton's  law)  for  the  system.  Since  weather  fluctuates  rapidly  we  determine 
an  effective  action  for  climate  by  integrating  over  all  weather  variations 

exp  - (effective  climate  action)  * Jd  (weather)  exp  - (total  action)  . 

This  leads  directly  to  the  effective  potential  point  of  view. 


A 


This  paper  has  repeatedly  fallen  back  on  the  ice  cap- insolation  energy 
budget  models  as  a guide  to  the  formulation  and  possible  arena  of  use  of 
the  effective  potential  formulation.  One  may  well  imagine  that  in  many 
cases  where  non-linear  phenomena,  stable  and  unstable,  are  driven  by 
external  forces,  such  as  insolation,  the  effective  potential  will  provide 
a convenient  and  tractable  strategy  for  the  study  of  their  time  dependence. 
The  birth,  growth,  and  decay  of  dust  storms  on  Mars  is  just  such  an  example; 
its  study  in  the  framework  outlined  here  will  be  given  in  another  publica- 
tion (Vesecky  and  Abarbanel,  1977). 
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